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ABSTRACT 

Optical nebular emission lines are commonly used to estimate the star formation rate of 
galaxies and the black hole accretion rate of their central active nucleus. The accuracy of the 
conversion from line strengths to physical properties depends upon the accuracy to which the 
lines can be corrected for dust attenuation. For studies of single galaxies with normal amounts 
of dust, most dust corrections result in the same derived properties within the errors. How- 
ever, for statistical studies of populations of galaxies, or for studies of galaxies with higher 
dust contents such as might be found in some classes of "transition" galaxies, significant un- 
certainty arises from the dust attenuation correction. In this paper we compare the strength 
of the predominantly unobscured mid-IR [Nell] A15.5/im+[Ne III] A12.8/im emission lines 
to the optical Ha emission lines in four samples of galaxies: (i) ordinary star forming galax- 
ies (80 galaxies), (ii) optically selected dusty galaxies (11), (iii) ULIRGs (6), (iv) Seyfert 2 
galaxies (20). We show that a single dust attenuation curve applied to all samples can correct 
the Ha luminosity for dust attenuation to a factor better than 2. Similarly, we compare [Oiv] 
and [Olll] luminosities to find that [Olll] can be corrected to a factor better than 3. This shows 
that the total dust attenuation suffered by the AGN narrow line region is not significantly 
different to that suffered by the starforming H II regions in the galaxy. We provide explicit 
dust attenuation corrections, together with errors, for [Oil], [Olll] and Ha. The best-fit aver- 
age attenuation curve is slightly greyer than the Milky- Way extinction law, indicating either 
that external galaxies have slightly different typical dust properties to the Milky Way, or that 
there is a significant contribution from scattering. Finally, we uncover an intriguing corre- 
lation between Silicate absorption and Balmer decrement, two measures of dust in galaxies 
which probe entirely different regimes in optical depth. 

Key words: Type keywords here. 



1 INTRODUCTION 

In the last few decades large optical spectroscopic surveys have 
dramatically improved our knowledge of the global distribution of 
galaxy physical properties such as stellar mass, star formation rates, 
star formation histories, and both gas phase and stellar metallicities. 
Furthermore, at low to intermediate redshifts, the sheer statistics of 
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current optical spectroscopic surveys allow us to identify and study 
significant numbers of rare objects which represent transient phases 
in the lifetime of galaxies: starburst, post-starburst and "green- 
valley" galaxies, major mergers, and powerful active galactic nuclei 
(AGN). These transient phases, although rare, are globally impor- 
tant. For example, bl ack hole growth occurs primarily in rare, hig h 
accretion rate events ( lYu&Tremaine1l2002l :l Heck man et al]|2004h . 
and it is still debated whether the infrequent gas-rich galaxy ma- 
jor mergers are responsible for the morphological transformation 
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of galaxies and the build - up of the red sequ ence dWild et alj|2009l ; 
lEliche-Moral et al.ll2010l : lllbert et alJIioiO) . In the near future, the 
James Webb Space Telescope (JWST) will allow rest-frame opti- 
cal spectroscopic observations at high-redshift, where galaxies are 
more active in both star formation and black hole growth. 

There is an obvious benefit in using the information provided 
in the redshift survey to study the properties of the galaxies in de- 
tail. For example, when star formation is the dominant mechanism 
of excitation, Ha is a very simple and reliable measure of star for- 
mation rate in galaxies. However, it must be accurately corrected 
for dust attenuation. In ordinary galaxies, the tight relations be- 
tween dust-corrected Ha and other indicators of star formation rate 
from other wavelength regimes provide circumstantial evidence 
that dust attenuation i n the optical can be corrected for accurately. 
iKennicutt et al. suggests that errors on SFRs caused by the 

assumed dust curve (attenuation as a function of wavelength) can 
be ±10-20% for typical star-forming galaxies in the lo cal Universe. 
Howev er, knowing the "correct" dust curve is crucial. lMuzzin et al.l 
(2009) find that the use of different dust curves in spectral-energy- 
distribution (SED) fitting of high redshift K-selected galaxies can 
alter parameters such as star formation rates and star formation 
timescales by factors of 2-3. The large errors reported in this study 
result from the fact that the "correct" dust-curve to use at high 
redshift is largely unconstrained and there is a large range in the 
steepness of dust-curves observed along lines-of-sight to stars in 
the Magellanic clouds and Milky Way. 

With the complex physical processes that occur during tran- 
sient phases in the life-cycle of a galaxy, it is uncertain how far 
rest-frame optical studies alone can help us in understanding their 
true impact. As dust content increases, and galaxy geometries be- 
come more complicated, corrections for dust attenuation become 
more uncertain. How much physical information can we really in- 
fer from optical data alone when faced with extreme star formation, 
dust production, AGN activity and complicated spatial geometries 
of disturbed galaxies? 

In this paper, we compare emission line fluxes in the opti- 
cal to those in the mid infra-red where dust attenuation is negli- 
gible, except in small and well known wavelength bands. To test 
the recovery of Ha luminosity from optical data in the presence 
of significant dust, we compare and contrast a low-z "representa- 
tive" sample of galaxies with (i) dusty galaxies selected from the 
Sloan Digital Sky Survey (SDSS) based on their Balmer emission 
line ratios and strong Balmer absorption lines indicating a possi- 
ble recent decrease in star formation; (ii) local Ultra-luminous In- 
frared Galaxies (ULIRGs) selected from the Infrared Astronomy 
Satellite (IRAS) Bright Galaxy Sample (BGS) based on their to- 
tal IR luminosities (Ltir > 10 12 ). The control sample of galaxies 
arise from the Spitzer-SDSS-GALEX Spectroscopic Survey (SS- 
GSS) and were chosen to be "representative" of the ordinary galax- 
ies within the SDSS and GALEX surveys. For the AGN in these 
samples, we further compare the [OlII]/[OlV] line ratio to a sample 
of Seyfert 2 galaxies, to test the recovery of black hole accretion 
rates from optical [OIII] emission line. In a companion paper we 
will address the effect of dust on the optical continuum, and thus 
estimations of star formation histories, and compare the star forma- 
tion rates derived from the multi-wavelength continuum with those 
from Ha luminosity. 

Understanding our selected "unusual" galaxies is crucial for 
understanding galaxy evolution. The B aimer-strong galaxies show 
optical emission lines indicating obscured AGN, and the correction 
of their [OIII] emission for dust attenuation using different dust- 
curves results in them contributing between an insignificant ~1% 



and a much more significant -~10% to the total [OIII] luminosity (or 
black hole accretion rate) of all narrow-line AGN in the SDSS. This 
has important implications for the starburst-AGN connection, for 
example, as the strong Balmer absorption lines sugge st this accre - 
tion is happening after the starburst has decayed (Wild et~alj2010h . 
Lo w-z ULIRGs are usu ally associated with major gas-rich merg- 
ers ^Sanders et alJI 1988b . Some models of galaxy evolution have 
used gas-rich major mergers as the driving force to change the ob- 
servational properties of the galaxy population over cosmic time. 
The mergers lead to powerful starbursts and QSO driven outflows, 
resulting in the po pulation of red, quiescent ellipticals s een in the 
local universe (e.g. lHopkins et alj2006l ; lNaab et alj2006h . ULIRGs 
become increasingly common at earlier epoc hs, and dominate th e 
star formation rate density beyond a z ~ 1 jBorgne et al.| [2009). 
Understanding the basic properties of ULIRGs at low-redshift, and 
the limitations of the optical wavelength regime for identifying and 
understanding the properties of their power sources, is crucial in an 
era when near-infrared (NIR) ground-based spectrographs and the 
IWST will soon allow us to study galaxies with rest-frame optical 
spectra at increasingly early epochs. 

The outline of this paper is as follows. After briefly reviewing 
the topic of dust extinction and attenuation in Section[2] we present 
the spectral analysis methods and four samples in detail in Section 
[3] In Section|4]we compare the data with a variety of different dust 
curves commonly used in the literature. We conclude in Section 
[5] by discussing the implications of our work for measuring star 
formation rates, black hole accretion rates, and identifying AGN 
in optical spectroscopic surveys. Finally, we provide relations for 
correcting the [Oil], [OIII] and Ha emission lines using a measured 
Ha/H/3 flux ratio (Balmer decrement), assuming a single dust curve 
which provides a good fit to all of the samples studied in this paper. 

Throughout the paper we quote line luminosities in erg/s and 
masses in solar masses. Optical parameters for all samples have 
been extracted from the SDSS-DR7, and may differ from those in 
earlier releases due to updates in the SDSS reduction pipeline. All 
optical wavelengths are quoted as measured in a vacuum, in keep- 
ing with SDSS convention, optical wavelengths are given in A and 
mid-IR wavelengths in fim. 



2 REVIEW OF DUST ATTENUATION LAWS 

Dust grains in galaxies which intercept and scatter stellar and nu- 
clear light complicate the analysis of integrated galaxy spectra. Op- 
tical light reaching the observer is effected by dust in three dif- 
ferent ways: (i) absorption; (ii) scattering out of the line-of-sight; 
(iii) scattering into the line-of-sight. The first two combined are 
usually termed "extinction" and all three combined are termed "at- 
tenuation". The final "effective attenuation" observed in the inte- 
grated light of galaxies results from the integration of attenuated 
light over all lines-of-sight to light sources throughout the entire 
galax)B In this section, we introduce the basics of the effects of 
dust on studying the integrated light from galaxies. ICalzettil J200lb 
gives an extensive recent review of the literature, although note that 
the terminology differs in places from that adopted here. 



1 Throughout this paper we use the term "effective attenuation" to encom- 
pass all effects relating to the dust i.e. absorption and scattering, and ge- 
ometrical effects. We note that this differs from some other papers which 
apply the term to "single-line-of-sight" attenuation. 
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The effective attenuation A\, in magnitudes, at a given wave- 
length A is given by: 



dust model: 



-2.5 log 



I\,ob 



= -2.5 log [exp(-n)] 



(1) 
(2) 



where I\, bs is the observed luminosity, Ix^ntr is the intrinsic lu- 
minosity of the source, and t\ is the effective optical depth of the 
dust. 

Observationally we measure the dimming and reddening of 
light from sources, thus deriving an effective attenuation or extinc- 
tion curve k\ : 

A\ 



kx = 



E(B - V) 
E(A - V) 
E(B - V) 



(3) 



(4) 



where the colour excess is defined by E(\ — V) = A\ — Ay and 
Rv is the total-to-selective extinction in the V band. In the case 
of extinction, Rv is a function of the grain composition, size and 
shape. In the case of effective attenuation observed in galaxies it is 
additionally a function of the local geometry of the dust and light 
sources and global geometry of the galaxy (bulge and disk, and 
inclination). In this work we will use the term "dust curve" to mean 
either "effective attenuation" or "extinction" curve, depending on 
the particular curve under discussion. 

The effective attenuation suffered by the ultra-violet (UV) to 
optical stellar continua of starburst galaxies, where the unattenu- 
ated stellar continuum is relatively well understood, has been mea- 
sured empiricall y from the ratio of heav ily attenuated to less atten- 
uated spectra JCalzetti et all 1994 [2000). The dust-curve is found 
to be "greyer" (i.e. flatter) than extinction curves measured along 
single lines-of-sight towards stars in the Milky W afl as expected 
from geometrical effects. Similarlv. lCharlot & FallH2000] . hereafter 
CF00) fit the ratio of far infra-red (FIR) to UV luminosities and UV 
slope in starburst galaxies to find an effective attenuation curve that 
is well fit by the shallow powerlaw t\ oc A -0,7 . 

Dust curves derived from stellar continua are not necessar- 
ily suitable for the correction of nebular emission lines, however. 
In starburst galaxies, nebular emission l ines are found to be more 
attenuated than the stellar continua (e.g. ICalzettietal]| 19941) . indi- 
cating that the light passes through a greater quantity, and possibly 
different distribution or composition, of dust. Any difference in the 
distribution or composition of dust will result in a different shape 
for the dust attenuation curve and therefore errors in the correc- 
tion of nebular emission lines using a dust-curve derived from the 
stellar continua. The stronger attenuation of nebular emission lines 
compared to the stellar continuum, can be successfully reproduced 
by models in which dust is distributed primarily in 2 components: 
diffuse interstellar dust, and stellar birthclouds which disperse after 
a finite time (~ 10 7 years) JCalzettietal]| 19941 ; Isilva et al.|[l998l. 
CF00). This allows the emission lines, which are ionised by the 
hottest stars still enveloped in their dense birthclouds, to experience 
greater attenuation than the UV stellar continua, which has an addi- 
tional contribution from cooler stars around which the birthclouds 
have already dispersed. 

In Wild et al. (2007) we presented explicitly the form of the 
dust-curve for the correction of nebular emission lines in the CF00 
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We reiterate that we use ry to refer to the final "effective atten- 
uation" suffered by light, including geometrical effects, and that 
this model is strictly an "angle averaged" model, averaging over 
all sight-lines through a galaxy. The first term describes the slope 
of the screen-like extinction suffered as the light leaves the stellar 
birth-clouds and the exponent is set to match the extinction ob- 
served along lines-of-sight in the MW (n = 1.3). The second term 
describes the additional, much greyer attenuation caused by the in- 
tervening diffuse ISM, where CF00 derive a best-fit exponent of 
m = 0.7. The fraction of the total effective optical depth at 5500A 
caused by dust in the diffuse ISM is controlled by the parameter fi. 
CF00 showed that the typical value for n in their starburst dataset 
was about 0.3. The resulting model attenuation curve, labelled /io.3 
throughout the paper, naturally lies between the pure continuum 
t\ oc A~ 0,7 attenuation (equivalent to /i = 1), and pure screen-like 
attenuation (e quivalent to fi = 0). 

Recentlv. lda Cunha et alj a2008h applied and extended the two 
component CF00 dust model to a sample of local galaxies with 
multi-wavelength photometry and optical Ha and H/3 emission line 
measurements, finding an excellent match to the data. 

We note that an additional complication arises in the case of 
nebular emission from the narrow line region (NLR) of an AGN, 
where a different dust geometry could lead to a different effec- 
tive attenuation curve than that found for nebular emission from 
Hll regions. Additionally, for galaxies with contributions to neb- 
ular emission from both star formation and an AGN NLR, the 
use of the Balmer decrement for the correction of NLR emission 
lines may be somewhat inappropriate. In general it is reasonably 
assumed that the Hydrogen recombination lines are dominated by 
the emission from H II regions in all but the strongest AGN. How- 
ever, this is not the case for higher ionization and higher excitation 
lines where emission from the NLR can dominate. Thus, while the 
Balmer decrement will show how much dust is in H II regions, there 
is no reason to believe that light from the NLR would pass through 
the same amount of dusQ. Therefore the extinction correction ap- 
plied to these other lines, such as [[OIII]] A5007A, may be incorrect. 
We investigate these effects in Section |4~4l 



2.1 Attenuation of nebular emission lines 

In order to determine the amount of attenuation suffered by op- 
tical nebular emission lines, ideally one would observe a pair of 
lines whose intrinsic ratio is independent of physical conditions and 
where one of the pair lies at a wavelength l ong enough that atten- 
uation is negligible (see also discussion in lOsterbrock & Ferlan j 
2006). Hydrogen recombination line ratios that include one line in 
the NIR wavelength range have been shown to be a powerful tool 
for constraining the dust attenuation curves for neb ular emission 
lines in galaxies dCalzetti et ai]|l996l : ICalzettill200ll) . The data are 
well fit by a screen-like MW extinction curve, which has since been 
favoured for the correction of emission lines for dust attenuation. 
These results suggest that a substantial fraction of the attenuation 
of emission lines must arise from dust in a foreground screen-like 
geometry (e.g. the birthclouds in the CF00 model). However, the 



2 Commonly used M W extinction cur ves include ISeatonl 
ll979hlCardelli et all ll989l) and lO'Donnel] dl994l) 



3 The difference in the intrinsic Balmer decrements between NLR and H II 
regions (due to different gas temperatures), will be small relative to the 
difference caused by attenuation. 
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samples studied to date have been small, containing only local star- 
burst galaxies. Extending this study to a larger sample of galaxies, 
including more ordinary starforming galaxies and dustier galaxies 
is the purpose of this paper. 

Sample sizes using ground based NIR observations are re- 
stricted primarily by the difficulty resulting from contamination 
of observations by OH sky lines. Hubble Space Telescope obser- 
vations of NIR Paschen-a line have been used to reveal the de- 
tailed distributio n of star format ion in galaxies relatively unim- 
peded by dust (Bok er et al. 1999). The detection of dust lanes in 
some galaxies in this sample indicate some level of extinction still 
exists even at 1.875^im. Uncertainties in the form of dust-curve in 
the NIR regime therefore cause a small additional source of uncer- 
tainty when using NIR emission lines to constrain dust-curves in 
the optical. 

With the launch of th e Spitzer Space Te lescope, with the In- 
frared Spectrograph (IRS, Houck et alj|2004h on board, high qual- 
ity, moderate resolution mid-IR spectra are available for the first 
time for a large number of low-z galaxies. Extinction in the mid-IR 
is negligible, except in small wavelength regions centered on 9.7 
and 18pm where strong absorption can arise from the stretching 
and bending modes of amorphous silicate grains. While H recom- 
bination lines are found in the Mid-IR, they unfortunately tend to 
be too weak to be easily measured. Two of the strongest and eas- 
ily observable emission lines in the mid-IR arise from neon; [Ne II] 
12.8/im and [Nelll] 15.5^m. In typical Hll regions, over 95% of 
neon is in these two states (Ne + and Ne 2+ ), thus we expect that 
the sum of these two lines wi ll trace very well th e ionized gas, and 
thus Ha, in galaxies. Indeed. iHo & Ketd have already sug- 

gested that these mid-IR Neon lines could provide a useful new 
indicator of star formation rate in galaxies. They have shown that 
[Ne Il]+[Ne III] luminosity correlates strongly with both total in- 
frared luminosit y in galaxies and the NI R Bra emission line in gi- 
ant H II regions. Diaz-Sant os et alj d2010t> have shown that [Ne II] 
luminosity correlates strongly with that of Paa in star-forming re- 
gions in LIRGS. Neon has the added benefit of being an abundant 
element, especially as it is not depleted onto dust grains. Therefore 
in this work we use the sum of these two lines as a tracer of the un- 
extinguished emission from ionized gas in these galaxies. Although 
these lines lie on the edge of the silicate absorption features, in gen- 
eral the extinction caused by this is minimal as we verify in Section 
[3] Combining Spitzer spectra with optical SDSS spectra provides 
a new and extremely powerful dataset, in which both optical and 
mid-IR emission lines have been measured much more accurately 
than was possible in the past. 

Starburst/H II r egion modelling dKewlev et alj 1200 ll : 
iLevesque et alj [2010) suggests that for typical twice Solar 
metallicity H II regions, the ratio L([Ne Il]+[Ne IIl])/L(Ha) should 
be between 0.5-1.0, demonstrating how luminous the mid-IR Ne 
lines are. The range is caused by ionization effects, arising from 
different stellar models, the mean age of the ionizing stellar popu- 
lation, and ionizing flux densities. There is one major disadvantage 
of using the [Ne]-to-Ha ratio to constrain the optical dust-curve 
in galaxies, in that the intrinsic ratio is directly dependent upon 
the metallicity of the gas, greatly expanding the range of possible 
values. For Solar metallicity Hll regions the models predict 
L([NeIl]+[NeIIl])/L(Ha) should be smaller by a factor of 0.6. 
However, we expect that this is not a dominant concern in the 
samples studied in this paper as we are exploring massive, dusty, 
and therefore metal-rich, galaxies. We discuss the impact of 
variations in the intrinsic ratio where relevant. We refer the reader 



to lHo&K eto (2007) for further details on the formation of the 
neon lines in H II regions. 

In AGN the ionizing spectrum is much harder than found in 
star-forming regions, meaning that neon starts to populate higher 
ionization states. The sum of the mid-IR neon lines still traces the 
ionized gas very well, but the ionisation is offset to higher values 
due to the harder radiation field. In addition, the harder radiation 
field of AGN also create partially (or X-ray) ionized regions, which 
have a higher L([[Ne II]])/L(Hq) ratio. AGN photoionization mod- 
els suggest L([N e Il]+[Ne IH])/L(Ho Q ~ 2, a factor of 2-4 higher 
than H II regions dGroves et alj2004TFl Again, we will indicate the 
effect of this result where relevant. 



3 THE SAMPLES 

In this section we give the methodology used to measure the 
nebular emission lines, followed by the selection criteria of our 
four samples of galaxies. All samples have SDSS optical 5- 
band imaging and fibre spectroscopy, together with Spitzer IRS 
spectroscopy. The SDSS specobjid, SDSS-MPA rowindex, RA 
and Dec of the B aimer-strong, ULIRG and Seyfert samples are 
given in Tables IA1I IA2I and I A3 1 to allow easy identification of 
the objects in the catalogues used in this paper. The equiva- 
lent information for the SSGSS sample is available online at 
http://www.astro.columbia.edu/ssgss For a quick visualisation of 
the samples: Figure Q]presents their distributions in redshift, stellar 
mass and Balmer decrement (Ha/H/3 flux ratio), Figure [2] shows 
their positions on t he BPT diagram ([N ll]/Hq vs. [OlII]/H/3 emis- 
sion line flux ratios Baldw in et alj 1981i), and in Figure[3]we place 
all our samples on the ISpoon et a lj J2007I) diagram of PAH3.2 
equivalent width vs. the strength of the 9.7/im silicate absorption. 
We discuss these figures further in Section |331 

The raw optical emission-line measurements for all SDSS- 



DR7 galaxies are provided on the SDSS- MP A web pagej fl and 
details are given in lKau ffmann et al. ( 2003) and Brinchman n et al.l 
j2004b : errors have been rescaled according to the information pro- 
vided on the web site. Note that all broad emission line objects 
have been identified by the SDSS automated pipeline and are not 
included within the SDSS-MPA sample. Briefly, spectral synthesis 
models [Bruzual & Chariot (2003); 2008 version using the MILES 
spectral library, Sanchez-Blazquez et al. (2006)] are used to fit the 
stellar continuum, including stellar absorption features, and emis- 
sion lines are measured from the residual of the data and model fit. 
Therefore underlying stellar absorption is accurately accounted for 
when measuring optical emission line fluxes. Intervening Galactic 
extinction is corrected for prior to the calculation of the emission 
line fluxes. The quality of the fits obtained is excellent. 

For all samples, IRS spectroscopy was carried out in staring 
mode and ID spectra were extracted using the standard Spitzer 
IRS Custom Extraction (SPICE). For the B a imer-strong galaxie s 
the "optimal extraction" option was chosen dNarron et al]|2007h . 
which weights the extraction by the object profile and the signal- 
to-noise of each pixel and is suitable for faint sources. Individual 
ID spectra were combined using custom algorithms. 

To maximise consistency between the heterogeneous sam- 
ples in this paper, the final low-resolution ID spectra of all sam- 
ples were reanalysed using the same custom algorithms and anal- 



4 These line ratio predicti ons can be obtained u sing the publicly available 
ITERA software package jGroves & Alien]|201(]|) 

5 http://www.mpa-garching.mpg.de/SDSS/DR7/ 



© 0000 RAS, MNRAS 000, 000-000 



Nebular emission from dusty galaxies 5 




Figure 1. The distributions of redshifts (left), stellar masses measured from SDSS 5-band optical photometry (center) and Ha/H/3 ratios (right) of the 4 galaxy 
samples: SSGSS (green, shaded), Balmer-strong (red, line-filled), ULIRG (purple), Seyfert (blue, dotted). The black line shows the distributions for SDSS- 
DR7 galaxies which have been spectroscopically observed and have a spectral SNR g > 8. The downward arrow in the central panel shows the position of the 
observed bimodality in the galaxy mass distribution (because the histograms are not corrected for sample incompleteness the bimodality is not noticeable in 
this figure). 



ysis packages where appropriate. To obtain [Ne II] A15.5/im and 
[Nelll] A12.8/im fluxes for the SSGSS and Seyfert samples we 
used the mid-IR spectral fitting package PAHfit (vl.2 lSmith et al.l 
l2007t) which fits Gaussian profiles with widths set according to the 
instrument module. PAHfit was run in the rest-frame using a screen- 
like dust model, and although dust attenuation is tiny at the wave- 
lengths of the Ne lines, for complete consistency betw een samples 
we re moved the correction applied by the code (Eqn. 4, Smi th et ail 
120071) . Our final results are unchanged whether the minimal dust 
correction is applied or not. Where high resolution spectroscopy is 
available, we use the high-resolution measurements of [Ne II] and 
[Ne III] due to the smaller error bars. For the Balmer-strong sample, 
the [Ne II] and [Ne III] emission line fluxes were measured interac- 
tively from the short-high (SH) IRS spectra using the 2-step base- 
line and Gaussian lin e-fit module from the SMART IDEA package 
(Higd on et al] [2004fl. For the ULIRG sample, fluxes were taken 
from h igh resolution spectroscopy, given in Table 4 of lArmus et al.l 
d2007h . To ensure consistency between the four samples, where 
possibl^l we verified that [Ne II] and [Ne III] fluxes measured with 
PAHfit or directly from the low-resolution spectra for the Balmer- 
strong and ULIRG samples were equal within the errors to those 
measured from the high-resolution spectra (with and without the 
minimal correction for attenuation applied in PAHfit). For the 
Seyfert sample, [Ne II] A15.5/im fluxes were found to match within 
the erro rs those measured in th e high-resolution mode, given in Ta- 
ble 6 of lLaMassa etai] fcOld) . 

[OIV] A26/im was measured from high-resolution (LH) IRS 
spectroscopy in all cases (measurement from the low-resolution 
spectra is not possible due to a nearby Fell line). For the Balmer- 
strong sample, [OIV] line fluxes were measured interactively using 
the SMART IDEA package as described above. For the ULIRG 
and Seyfert samples [OI V] line fluxes were ta ken from the tables in 
Annus et al. (2007) an d lLaMassa etallfeOlti) respectively. 

One of the major components of a star-forming galaxy spec- 
trum in the mid-IR is the broad emission features generally ac- 
cepted to result from the vibrational modes of polycyclic aromatic 
hydrocarbons (PAHs). While we do not study the PAH features 
in detail in this paper, PAHfit outputs their fluxes and equivalent 



6 http://isc.astro.cornell.edu/SmartDoc/LineFitting 

7 We found that PAHfit does not fit the continuum well for ULIRGs with 
strong Silicate absorption, and therefore the emission line measurements 
were unreliable in these cases. 



widths which we use to provide further clarification of the nature 
of our objects where necessary. 

To visualise the different dust and AGN properties of our dif- 
ferent samples in Section [331 we measure the equivalent width of 
PAHs. 2 and the strength of the 9.7/im silicate absorption. We at- 
tempted in both ca ses to follow closely the methods as described in 
ISpoon et all ( 120071) . and we used the ULIRG s ample to ver i fy tha t 
our results are consistent with those given in ISpoon et al.l bOOTh 
over the full range of the measurements. 

3.1 The SSGSS sample 

Our control samp le is the Spitzer-SDSS-GALEX Spectrosco pic 
Survey (SSGSS, lO'Dowd et al] 120091 : Ijohnson et all |2009|) , a 
Spitzer Legacy survey (P.I. D. Schiminovich) which targeted 100 
"ordinary" galaxies from the SDSS spectroscopic survey and 
GALEX surveys for follo w-up mid-IR spec troscopy with the In- 
frared Spectrograph (IRS. [Houck et al] |2004|) using the Spitzer 
Space Telescope. The galaxies are magnitude limited to the SDSS 
spectroscopic magnitude limit of r < 17.7 and additionally were 
selected to have GALEX-NUV fluxes > 1.5mJy for the faint sam- 
ple (40 galaxies) and > 5mJy for the bright sample (60 galaxies). 
They were chosen such that their mass, colour, star formation rate 
and redshift distribution are representative of galaxies in the SDSS 
survey. A description of the survey, together with data tables con- 
taining SDSS and mid-IR measurements and derived parameters, 
are available at the team websitjf). In this paper we make use of 
the IRS low resolution spectroscopy of 80 SSGSS galaxies. The 
IRS spectra of the remaining galaxies are either incomplete due to 
problems during observations or reductions, or suffer from consid- 
erable noise at short wavelengths, 5 of these caused by a solar flare. 
A further 4 galaxies with median per-pixel-SNR< 10 in their SDSS 
spectra are removed, as the spectra were of insufficient SNR to al- 
low accurate stellar continuum fitting around the B almer emission 
lines. Furt her details of the SSG SS can be found in lO'Dowd et al] 
d2009h and ljohnson et al] ]2009l) . 

3.2 The dusty Balmer-strong sample 

Our first sample of galaxies with significant dust contents are 12 
Balmer-strong galaxies selected from the SDSS-DR4 to have an 
observed Ha/H/3 flux ratios greater than 8 (where a ratio close 

8 http://www.astro.columbia.edu/ssgss/ 
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to 3 indicates little dust) and optical line ratios that indicated the 
possible presence of an AGN. Additionally, they were selected to 
have strong Balmer absorption lines, which might indicate a recent 
sharp decline in star formation rate ("post-starburst"). However, the 
strong Balmer absorption lines might alternatively be caused by 
unusual dust properties rather than star formation history(e.g. as 
postulated in IPoggianti & Wulh oOO). This is the subject of a future 
paper. Both high- and low-resolution IRS spectra were obtained for 
all of these 12 galaxies during Spitzer Cycle-4 (P.I. T. Heckman). 
Further observational details of this sample are provided in Ap- 
pendix [B] One of the 12 galaxies was subsequently discovered to 
have double peaked narrow Ha and [Nil] lines, together with sig- 
nificant underlying broad emission which was not accounted for in 
the original emission line measurements. Its IRS spectrum is dom- 
inated by emission from the AGN. Attempts to robustly remeasure 
the H/3 line accounting for the double peak were unsuccessful due 
to significant stellar absorption, and its Ne/Ha and PAH/Ha ratios 
suggests that it is not in fact unusually dusty. We therefore remove 
it from all further discussion, leaving 1 1 Balmer-strong galaxies. 

3.3 The ULIRGS 

Our second sample of dusty galaxies is composed of the 10 
ULIRGs (Lir > 10" L B ) from the IRAS Bright Galaxy Sam- 
ple jSoiferetalj [T987h: the Spitzer IRS low- and high-resolution 
ob servations of these galaxies have been presented in detail 
in lArmus et alj fc007l) . Of the 10 BGS ULIRGs, 6 are found 
in the SDSS spectroscopic sample [IRAS 15327+2340 (Arp220), 
IRAS 08572+3915, IRAS 15250+3609JRAS 09320+6134 (UGC 
5101),IRAS 12112+0305], 3 do not lie within the SDSS footprint 
(IRAS 05189-2524, IRAS 22491-1808, IRAS 14348-1445), and 1 
is optically too faint to have been targeted [IRAS 12540+5708 (Mrk 
231) has an SDSS r-band magnitude of 1 9 4], All 6 of the ULI RGs 
with SDSS spectra are late stage mergers dVeilleux et alj2002h . 

3.4 The Seyfert sample 

In Section l4!4l we compare the high ionisation emission lines in the 
dusty galaxies with those of a sample of relatively dust-free (low 
Balmer decrement) Seyfert galaxies. This is composed of the top 20 
[OIII] A5007 flux emitters out of all SDSS-DR4 main-sample galax- 
ies classified as Seyfert galaxies based on their position on the BPT 
diagram. They have [OIII] fluxes greater than 4 x 10~ 14 erg/s/cm 2 , 
corresponding to a luminosity of ~ 10 41 erg/s at the median redshift 
of the sample. Both high- and low-resolution IRS spectra were ob- 
tained, with high signal-to-noise ratio. Further details of the Seyfert 
sa mple, along wi t h a stu dy of their X-ray proper ties can be found 
in lLaMassa et al.1 d2009t) and lLaMassa et alj ( bold) . 

3.5 Sample comparison 

The samples that we will compare and contrast in this paper were 
targeted for different purposes and thus any conclusions drawn 
must take account of potential systematics caused by their differ- 
ent selections. 

The redshift distribution of the 4 samples is shown in the left 
panel of Figure[T] The black line shows the redshift distribution of 
the SDSS spectroscopic galaxy catalog with r-band extinction cor- 
rected petrosian mag nitudtJE] in the range 14.5 < r < 17.77. The 

9 Petrosian magnitudes and radii are estimates of total mag- 



log N 




log f[NII]/f(Ho) 

Figure 2. The BPT diagram I Bald win et aT] Il98ll) : the narrow 
emission line flux ratios BPT^ =f[NII],\6585/f(H a )A6565 and 
BPT,; =f[Offl] A5 oo8/f(H/3)A4863 ft* the SDSS-DR7 spectroscopic 
galaxy sample (gray scale), SSGSS starforming galaxies (green dots, 
35/35), SSGSS AGN (green open triangles, 4/4), SSGSS composite 
galaxies (green open diamonds, 21/21), Balmer-strong galaxies (red dots, 
11/11), ULIRGs (purple stars, 6/6), Seyferts (blue triangles, 20/20). Only 
galaxies with all 4 lines detected at > 3cr confidence are plotted. The 
dashed lines indicate our division between starforming/composite (blue) 
and composite/AGN (black) samples. These demarcation lines define the 
three optical classes (starforming, composite, AGN). To aid the discussion, 
the Balmer-strong galaxies and ULIRGs in the AGN region have been 
labelled. 

Seyfert and ULIRG sample are flux limited in nebular emission and 
IR flux respectively, and lie predominantly at low redshifts. The SS- 
GSS galaxies were selected to be representative of the full SDSS 
sample, and therefore their redshift distribution tracks closely that 
of the full SDSS survey. The Balmer-strong sample is limited to 
z < 0.07 to ensure that the stellar population sampled by the opti- 
cal fibre through which the spectrum is taken lies within the central 
region of the galaxy close to the AGN. 

The stellar mass distribution of the 4 samples is shown in 
the central panel of Figure Q] The stellar massesj 1 ^! have been esti- 
mated from stellar population model fits to the optical 5-band SDSS 
p hotometry, using a Bayesian a nalysis similar to that described 
in lKauffmannetal.1 j2003l) and lSalim et alj d2005h . Most samples 
cover the full range in stellar mass probed by the SDSS, although 
the Balmer-strong galaxies do not cover the low-mass end of the 
distribution. Note that the histograms are not corrected for incom- 
pleteness effects due to sample or survey flux limits, and therefore 
are not suitable for drawing conclusions about the underlying distri- 
butions of galaxy mass. To aid the reader, the position of the galaxy 
mass bimodality is marked by a downward arrow. The right hand 
panel of Figure [T] shows the Balmer decrement distributions of the 
4 samples, showing how the Balmer-strong and ULIRG samples 
are significantly dustier on average than the SSGSS and Seyfert 
samples. 

In Figure [2] we show the position of the four samples on the 

nitudes suitable for all types of galaxies. They are derived 

by the SDSS photometry pipelines, with details found here: 

http://www.sdss.org/dr7/algorithms/photometry.html 

lu Courtesy of Jarle Brinchmann, publicly available at 

http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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Figure 3. The strength of the 9.7/xm silicate absorption vs. the equivalent 
width of the 6.2/jm PAH fea ture, following the measurement methods pre- 
sented in ISpoonetalJj2007l) . Note that the x-axis is linear in order to better 
distinguish our samples. The dotted lines are two of the dividing lines given 
by Spoon et al. (2007) to separate galaxies with no significant Silicate ab- 
sorption (horizontal) and ordinary star-forming/starburst galaxies (vertical). 
Median errors are shown in the bottom left for the Seyfert and SSGSS sam- 
ples, median errors are smaller than the symbol size for the Balmer-strong 
and ULIRG samples. Only those galaxies with mean per-pixel-SNR> 2 be- 
tween 5.4 and 5.9/im are plotted. Symbols are: SSGSS starforming galax- 
ies (green dots, 25/35), SSGSS AGN (green open triangles, 2/4), SSGSS 
composite galaxies (green open diamonds, 14/21), SSGSS galaxies with in- 
sufficient SNR in one or more optical emission lines to provide an optical 
BPT classification (green crosses, 20/20), Balmer-strong galaxies (red dots, 
11/11, labelled), ULIRGs (purple stars, 6/6, labelled), Seyferts (blue trian- 
gles, and upper limits marked by blue arrows, 13/20). 



BPT diagram, with the underlying grayscale showing the distri- 
bution of spectroscopically observed galaxies in the SDSS. Only 
those galaxies with SNR> 3 in all four emission lines are plotted. 
This is the most commonly used optical diagram for diagnosing 
the presence of an AGN. Throughout this paper, we distinguish 4 
classes of galaxies: (i) AGN -dominated, with emission-line ratios 
w hich place them abov e the stringent theoretical starburst criterion 
of lKewlev etailfeOOlh . or have [Nll]/Ha> 0.0; (ii) star-formation 
dominated, w hich lie below the observ ationally determined demar- 
cation line of iKauffmann et all {2003), or have [NII]/Hq< -0.5; 
(iii) composite objects, which lie between the AGN-dominated and 
star-forming dominated samples; we term them "composite" to in- 
dicate that the emission may be caused by a co mbination of both 
AGN and starformation ( Kauffmann e t alj|2003l although see cau- 
tion below); and (iv) unclassified objects with non-existent or weak 
(SNR< 3a) emission lines which do not allow classification into 
the first three samples. 

We note that the BPT diagram does not provide a firm de- 
marcation between the samples and that the demarcation lines are 
not motivated by robust theoretical predictions; we use the lines 
only to provide a convenient broad classification. iLevesque et al] 



present up-to-date modelling of the "maximum starburst" 
demarcation line, placing it muc h closer to the obser ved star form- 
ing branch in the SDSS than the lKewlev et al] l l200ll) line. Shocks, 
high metallicity, and ionisation by old stars can also ca use galax- 
ies to lie in the "composite" class (Stasiriska et al. 2008). We note 
that most of the Balmer-strong galaxies have detectable mo lecular 
H 2 lines, as do the majority of ULIRGs faigdon et alj200fj|) . In the 
case of ULIRGs these lines are thought to arise primarily from pho- 
todissociation regions (PDR) associated with star formation. The 
H2/PAH flux ratio of the Balmer-strong galaxies is exactly that ex- 
pected for star forming galaxies (Ogle et al. in prep.), indicating 
that shocks are similarly not a dominant excitation mechanism in 
the Balmer-strong galaxies. This suggests that the high [Nll]/Ha 
line ratios in these samples are in general caused by a central AGN, 
rather than shocks. Signal-to-noise ratio limits placed on emission 
lines may preferentially cause very high-metallicity galaxies and 
LINERs caused by shocks and/or ionisation by old star s to be "un- 
classified" in the SSGSS sample iFernandes et al. 2010). However, 
none of these caveats affect the results of this paper which do not 
rely on a precise identification of weak AGN, or on the precise po- 
sitioning of pure AGN or starforming galaxies on the BPT diagram. 
Where possible we include unclassified SSGSS galaxies in the fig- 
ures. 

Out of the 80 SSGSS galaxies, 4 are AGN-dominated, 21 
are composite objects, 35 are star-forming and 20 are unclassified. 
Out of the 11 Balmer-strong galaxies 2 are AGN-dominated and 
the remainder composite objects. 3 ULIRGs are AGN-dominated 
(Arp220, Mrk273 and UGC5101) and 3 are composite objects 
(08572, 12112, 15250). All Seyferts clearly lie in the AGN dom- 
inated region by definition. Figure [2] supports our assumption that 
the intrinsic L([NeIl]+[NeIIl])/L(Ha) ratio will not vary greatly 
between our samples. Most galaxies in the samples lie in the 
solar and greater region of this line r atio diagnostic diagram 
dGroves et al]|2006l ; ILevesque et al]|2010l), 

Finally, in Figure [3] we follow ISpoon et al] bOOTl) and plot 
the equivalent width (EQW) of the 6.2/im PAH feature against the 
strength of the 9.7/im silicate absorption. We plot only those ob- 
jects with mean SNR-per-pixel> 2 in the wavelength regions used 
to define the continuum for the measurement of the silicate absorp- 
tion at 9.7/im and with total flux/error> 3 in the 6.2/im PAH fea- 
ture. The lower branch is interpreted as a mixing between AGN and 
starformation activity, and the upper branch as an evolution from a 
dust obscured, to dust unobscured nuclear starburst or AGN (see 
Rowan-Robinson & Efstathiou 2009). Because of its focus on both 
dust and AGN activity, Fig.[3]is particularly useful for characteris- 
ing our 4 different samples: 

Galaxies in the SSGSS occupy the bottom right corner of the di- 
agram as expected for "ordinary" star-forming galaxies or weak 
narrow-line AGN. They have strong EQW(PAHe.2) and weak sili- 
cate absorption. The composite galaxies in this sample have a lower 
mean EQW(PAHb.2) than the star-forming galaxies, strengthening 
the argument that these galaxies contain an AGN. The large scatter, 
and significant number of objects with apparent silicate emission, is 
indicative of the generally lower quality of these spectra compared 
to our other samples. SSGSS galaxies with optical emission lines 
which are too weak to classify them on the optical BPT diagram do 
not differ in mid-IR properties from their classified counterparts. 

The Balmer-strong galaxies generally have stronger Silicate ab- 
sorption than SSGSS galaxies, placing them between the ULIRG 
and SSGSS samples and therefore, within the common interpre- 
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taion of the diagram, at the end of their evolution from a dusty star- 
burst. M82 is the prototype galaxy for this region of the diagram. 

The ULIRGS are a mixed group, as has been discussed at length 
in the literature. Although our sample is small, it is interesting 
to note that the ULIRGs with particularly strong silicate dust ab- 
sorption (08572, 15250), show optical emission line ratios indica- 
tive of a mixture of star formation and AGN and/or shocks. Three 
out of the four ULIRGs with weaker silicate absorption (Arp220, 
Mrk273, UGC5101), have AGN-dominated optical emission (i.e. 
implying little contribution from star formation, except possibly in 
the case of extreme shocks). This is consistent with the interpreta- 
tion of the upper branch as an evolution from dusty starburst galax- 
ies in the top-left to galaxies undergoing less intense, less dusty star 
formation in the bottom-right, with dust-obscured AGN activity oc- 
curring throughout the evolution. 

The Seyferts are generally found to have very weak PAH equiv- 
alent widths and little Silicate absorption, and therefore lie in 
the region of the diagram commonly interpreted as being oc- 
cupied by strong AGN which heat the dust to high tempera- 
tures. 2/20 show some silicate absorption together with stronger 
EQW(PAH 6 . 2 ). The median (mean) BPT„ values of the 4 Seyferts 
with EQW(PAH 6 . 2 )> 0.2 are 0.78 (0.75), compared to 1.01 (0.99) 
for the remainder of the sample. This suggests that emission from 
star formation contributes to the optical emission lines even in the 
Seyfert branch of the BPT diagram. These results are consistent 
with the interpretation of the lower branch of Fig. [3] as mixing be- 
tween AGN and starformation activity. 

3.6 Aperture Corrections 

'Aperture bias" describes the effects caused by the fact that a fixed 
aperture or slit width may probe only a small fraction of a galaxy's 
total light, and this fraction is a function of galaxy radius and red- 
shift. The physical size probed by the SDSS 3" fibre varies from 
under lkpc for the majority of the ULIRG sample to around lOkpc 
at the maximum redshift of the SSGSS galaxies. Fortunately the re- 
sults in this paper rely upon line ratios or equivalent widths rather 
than absolute luminosities, and therefore aperture bias in this usual 
sense is less relevant. 

However, comparing line strengths from the IRS spectra with 
those from the SDSS spectra leads to other aperture related issues. 
The finite (and different) slit and fibre apertures may mean that 
the IRS and SDSS spectra sample a different fraction of the same 
galaxy, and this will cause measured line ratios to differ from the 
true values. There are two important issues: (i) light is lost from 
the slit/fibre due to the PSF; (ii) light is lost from the outer regions 
of extended sources, depending on the size of the aperture and the 
accurate positioning of the aperture on the source. In this paper the 
precise flux calibration between the two datasets is relatively unim- 
portant, but it is still important that we minimise any additional 
scatter in the observed relations. To help the reader to understand 
the limitations of our data, we outline here the flux calibration pro- 
cedures of the mid-IR and optical spectra. 

The IRS PSF increases substantially from the blue to the red 
wavelength extremes of the instrument, and therefore each IRS slit 
has a different width de sig ned so as not to loose too much of the 
light at the reddest end3 lx l The spectral extraction procedure ac- 

11 The short-low (SL) IRS module covering the wavelength range 5.2- 
14.5^m has a slit aperture of 3.6" x57". The long-low (LL) IRS mod- 
ule covering the wavelength range 14.0-38.0/^m has a slit aperture of 



counts for the varying PSF along the slit by normalising the flux at 
each wavelength according to the fraction of light missed from a 
calibration stellar point-source. 

The SDSS DR7 spectra have been spectrophotometrically cal- 
ibrated using the PS F magnitudes of stars observed on the same 
plate as the galaxies jAdelman-McCarfhv et all2 008 ) . To aid aper- 
ture corrections, the SDSS-MPA emission line measurements have 
been recalibrated to the 3" fibre magnitudes. Therefore, to match 
the stellar PSF-normalised mid-IR spectra, we have removed this 
recalibration, returning the flux level of the SDSS emission lines to 
the stellar PSF-normalised values. 

The use of stellar PSF spectrophotometric calibration for both 
optical and mid-IR spectra allows good relative calibration for point 
sources, but inaccuracies will still arise in the case of extended 
sources. For centrally concentrated sources, such as circumnuclear 
starbursts in ULIRGs or AGN NLR emission, the point-source ap- 
proximation is a good one. However, for extended sources, the mea- 
sured mid-IR to optical line ratios may depend on the relative size 
of the apertures, and the relative surface brightness distribution of 
the sources. Both effects are difficult to correct for accurately, espe- 
cially in the case of radial differences within the sources (for exam- 
ple, in AGN vs. starforming nebular emission). We have therefore 
not attempted any further aperture correction, and will caution the 
reader at points in the text where this effect may be relevant. 



4 RESULTS 

In Figure |4] we plot Ha vs. [Ne II] A15.5^m+[Ne III] A12.8/im 
(from now on [Ne]) luminosities for all of the samples, firstly with 
Ha uncorrected for dust attenuation, and then corrected with five 
different dust-curves selected from the literature to cover the full 
range in possible greyness. The panels are ordered by increasing 
greyness of the dust-curve used, from the screen-like extinction- 
curves (SMC and MW), to galaxy continuum derived attenuation- 
curves (Calzetti and A" 7 ). The SMC curve was evaluated from 
the tabulated results o f IPeil dl992h with an R v = 3.1; the MW 
curve was taken from O'Donnel l (j 19941) with an Rv = 3.1; the 
Mo. 3 curve is from IWild et alj J2007I) ; the Calzetti curve is from 
ICalzetti et all feOOd) with R v = 4.05; the A" ' 7 curve is from 
CF00 as discussed in Section[2] 

In each panel a maximum-likelihood best-fit linear relation is 
fit to the SSGSS-starforming sample allowing for errors in both x 
and y axis quantities and fixing the intercept to be zero (dashed 
line). Twice and half the best-fit [Ne] luminosity for each Ha lu- 
minosity is indicated as dotted lines to aid comparison between 
the panels of the scatter between the two quantities. The best-fit 
gradient (g, equal to the best-fit intrinsic [Ne]/Ha ratio) and Pear- 
son's correlation coefficient (r) are given in the bottom right of each 
panel, both evaluated using the SSGSS-starforming sample only. 

In panel a it is immediately evident that without correcting Ha 
for dust attenuation the scatter between Ha and [Ne] emission lu- 
minosities is large, even for the SSGSS galaxies with normal levels 
of dust. For the dusty Balmer-strong and ULIRG samples, there is 
a clear excess of [Ne] for their observed Ha luminosity as expected 
for dusty galaxies. 

As we move through the dust laws, from least to most grey 
[panels b to/], we see steady trends. The best-fit intrinsic [Ne]/Ha 

10.6" X 168". At high resolution, the short-high (SH) module covers 9.9- 
19.6/im with a slit aperture of 4.7" X 1 1 .3" and the long-high (LH) module 
covers 1 8.7-37.2 with a slit aperture of 1 1 . 1 " X 22.3" . 
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Figure 4. The optical Ha>,g5g5 vs. mid-IR [Ne II] A15.5/im+[Ne III] A12.8^tm luminosities in erg/s. Pane/ (a): with no correction of Ha for dust attenuation 
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(green open diamonds, 19/21), Balmer-strong galaxies (red dots, 11/11), ULIRGs (purple stars, 6/6), Seyferts (blue triangles, 19/20). 



ratio monotonically decreases, and the correlation coefficient is 
maximal for the MW, /io.3 and Calzetti dust-curves. Focussing 
on the dusty Balmer-strong and ULIRG samples, we notice that 
these move from lying above the relation defined by the SSGSS- 
starforming galaxies in panels a and b, to lying below the relation 
in panel /. This is caused by their under- and over-correction for 
dust attenuation using the different dust-curves. 

In Figure[5]we show the samples more clearly by plotting the 
residual of the [Ne] luminosity from the linear fit using the MW, 
/io.3 and Calzetti dust curves. The mean and variance of each sam- 
ple in these figures is informative for comparing the different dust 
curves. For the dusty ULIRGs the mean is 2.0, 1.6, 1.5 for the MW, 
pio.3 and Calzetti dust-curves respectively i.e. they are significantly 
offset for all dust-curves (see below for further discussion). For the 
Balmer-strong sample, the mean is 1.1, 0.9 and 0.8, suggesting per- 
haps that a dust curve in between the MW and /10.3 is most appro- 
priate but any offset between this extreme dusty sample and ordi- 
nary starforming galaxies is very small. For the SSGSS starforming 
sample, the mean is 1 by definition for all dust-curves. The variance 
is ^0.15 for all dust-curves, i.e. the scatter is less than a factor of 2 



and there is no strong trend with Ha luminosity for the starforma- 
tion dominated galaxies. 

From these figures we can identify the expected offset of the 
AGN caused by a higher intrinsic [Ne]/Ha ratio. For the Seyfert 
sample and dust curve the mean offset is 1.26, i.e. about a 
25% enhancement in [Ne] for their Ha luminosity. The magnitude 
of this offset does not correlate with position on the BPT diagram or 
PAH equivalent width, although the size of the sample is too small 
to rule out the possibility that those with a smaller offset have a 
greater contribution to star formation in their lines. Averaged over 
the samples, the offset is seen only for the Seyfert galaxies, and the 
three SSGSS-AGN which lie in or close to the Seyfert region on the 
BPT diagram. The SSGSS composite sample have a mean of 1 for 
all dust-curves, which suggests that the AGN in general does not 
contribute significantly to the Ne/Ha line ratio in these composite 
galaxies. However, the trend upwards at low Ha luminosities is 
driven by composite objects, suggesting there may be a small effect 
at low star formation rates. The observed average enhancement in 
the [Ne]/Ha ratio of the ULIRGs noted above, could be due to 
small number statistics, the dust curve (i.e. a shallower dust curve 
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Figure 5. The residual [Ne] luminosity from the best-fit linear relation between dust corrected Ha and [Ne] luminosities, as a function of Ha luminosity in 
erg/s. From left to right, Ha has been corrected for dust attenuation using the MW, fXQ.3 and Calzetti dust-curves. The offset of the AGN-dominated sources 
from the starformation-dominated sources is clear. The mean offset for the Seyfert sample is a factor of 1.26 for the dust curve. Symbols and lines as in 

Fig.m 



would be more appropriate), or the presence of deeply buried AGN. 
We note that the ULIRGs are only offset by at most a factor of 3, 
i.e. dust corrected Ha luminosity tracks mid-IR [Ne] luminosity to 
better than a factor of 3 even in extreme dusty and disturbed objects. 
We will return to these points below. 

The maximum likelihood fits provide a measure of the dust- 
free ratio of [Ne] to Ha luminosity in our IRS/SDSS datasets, when 
star formation is the dominant excitation mechanism for the line 
emission: 

L[Ne] = g x L(Ha) corr (6) 

where L[Ne] = L([Nell] A15.5^m)+L([Ne III] A12.8^m) and g = 
0.89, 0.78, 0.75 for the MW, ^0.3 and Calzetti dust-curves respec- 
tively. While it is encouraging that the measured intrinsic ratio lies 
exactly at the center of the range predicted by the models for twice 
Solar metallicity Hll regions (see Section l2~T! . a precise compari- 
son between models and data would only be possible after much 
greater care has been taken to perform precise aperture correc- 
tions. In the case of the AGN, the emission line source is point-like 
and therefore the relative calibration between the mid-IR and opti- 
cal spectra should be the most accurate. The fact that the average 
[Ne]/Ha ratio is a little lower than predicted may be due to con- 
tamination by star formation in a subset of the Seyfert sample (see 
Sectionl331i. 

In the following subsection we study the differences between 
the dust curves in more detail. 



4.1 Comparison of common dust attenuation laws 

In Fig.[6]we plot the observed [Ne]/Ha vs. Ha/H/3 luminosity ra- 
tios for all of our samples. In the left-hand panel we use a loga- 
rithmic scale to allow the dustiest ULIRGs to appear. In the right- 
hand panel we focus on the less extreme objects, and only include 
galaxies where their optical emission line ratios indicate starforma- 
tion is a significant emission line source (classified as composite 
or starforming in the optical BPT diagram, Fig. [2}. Overplotted 
as black lines are the predicted line ratios from the different dust 
curves presented in the previous subsection. In each case the in- 



trinsic [Ne]/Ha ratio is set from the fit to the SSGSS-starforming 
galaxies as described above. 

Although there are several factors which can move galaxies 
within this diagram, dust has the strongest effect, increasing both 
the Balmer decrement and observed [Ne]/Ha ratio by an order of 
magnitude over the range probed by our samples. The different 
dust-curves predict different relations between these two line ra- 
tios, even when we allow the intrinsic [Ne]/Ha ratio to be a free pa- 
rameter fit by the data. As shown in the previous section, the MW, 
jito.3 and Calzetti dust-curves provide the best fits to the SSGSS- 
starforming galaxies. Figure[6]shows that the H0.3 and Calzetti at- 
tenuation curves also pass through the center of the B aimer-strong 
sample. 

The scatter of the SSGSS-starforming galaxies about the best- 
fit dust-curves is large, and this may have several origins. Firstly, 
the errors on the line flux measurements, particularly in H/3 are sig- 
nificant. Secondly, we expect some variation due to different metal- 
licities. The grey shaded region shows the range predicted by the 
theoretical models for twice Solar H II regions, and the arrow in the 
left hand panel shows how this region moves for Solar H II regions. 
The 3 galaxies with the lowest measured [Ne]/Ha ratio lie signif- 
icantly to the left of the 2Zq prediction, but within the Zq pre- 
diction. These 3 galaxies also have some of the highest [OlII]/H/3 
and [Ne IIl]/[Ne II] ratios of the sample, indicative of harder ioni- 
sation perhaps resulting from lower metallicity star formation. The 
single Seyfert galaxy with a [Ne]/Ha ratio less than 1 has the low- 
est [Nll]/Ha ratio of the whole sample, unusually low for Seyfert 
galaxies, and likel y associated with accretion of low-metallicity gas 
jGroves et al]|2006l) . 

While low metallicities can explain the scatter of SSGSS- 
starforming galaxies to the left of the dust-curve predictions, the 
scatter to the right is unlikely to be caused by metallicity variations 
as scatter in this direction requires significant and unlikely neon 
abundance increases (i.e. greater than twice solar). This leads us 
to the third, and most interesting, cause for scatter in the line ratio 
relation: intrinsic differences in the type or distribution of dust in 
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Figure 6. The observed [Ne]/H« and Ha/H/3 flux ratios. Left: A logarithmic scale allows all samples to appear, and all starforming, AGN, composite and 
unclassified objects are plotted. Right: A linear version of the same plot, with the axes restricted to focus on the less extreme galaxies. Only star-forming and 
composite galaxies are plotted. Overplotted as blac k lines are the predicted ratios from different dust curves: the SMC extinction c urve with an Ry = 3.1 
IPeil 19921) : a MW extinction curve with Ry =3.1 jO'D onnell 1993); me MO. 3 attenuation curve for emission lines I Wild et al. 2007); the Calzetti attenuation 
curve for starburst continua with Ry = 4.05 ICalzetti et alj2000l) ; t he A~° 7 attenuation curve for st arburst continua from CF00. The grey lines additionally 
mark a different MW extinction curve with Ry = 3.2 (dashed grev. lSea ton 1979; Nandv et all 19751) and our new best-fit dust-curve from Section l4~2l (solid 
grey) with fi = 0.4. In all cases the intrinsic Ne/Ha ratio is taken from the best-fit to the SSGSS-starforming galaxies. The grey shaded region shows the 
theoretical range in observed [Ne]/Ha ratio for a twice solar stellar population and the /xo.3 dust-curve, and the bottom arrow in the left hand panel indicates 
the theoretical decrease in this ratio as metallicity decreases from twice Solar to Solar. The upper arrow in the left hand panel indicates the average observed 
increase in [Ne]/H« ratio where an AGN is the dominant ionisation source (i.e. a factor of 1.26 observed for the Seyfert galaxies). For clarity, error bars are 
indicated for the ULIRGs in the left panel, and for the Balmer-strong galaxies in the right panel. Median errors are indicated in the right panel for the SSGSS 
starforming and composite galaxies. Symbols as in Fig. [4] 



galaxiej^l. In a 2-component dust model, increasing the fraction of 
dust in the diffuse ISM relative to the dense birthclouds increases 
the contribution from scattering and makes the dust-curve greyer. 
In the CF00 emission line dust model given in Eqn. [5] this rela- 
tive amount of diffuse vs. dense dust is controlled by the parameter 
/j,. By increasing fi to 0.9 (i.e. 90% of the effective optical depth 
from the dust comes from the ISM, and 10% from the stellar birth- 
clouds), we can describe the lower envelope of the SSGSS galaxies. 
While this qualitative result is suggestive, a further detailed study 
including orientation, metallicity and ionisation effects would be 
required to confirm that these galaxies do indeed have a signifi- 
cantly different dust geometry. 

It is evident that the 3 ULIRGs classified as composite galax- 
ies by their optical line ratios also appear to favour a greyer dust 
curve than the dusty Balmer-strong galaxies. However, as we dis- 
cuss further below, it is also plausible that the high Ne/Ha ratio 
observed in these ULIRGs, as well as in some SSGSS-composite 
galaxies, is simply due to a deeply buried AGN. 

We summarise the results up to this point by noting that a sin- 
gle appropriate dust curve, combined with an accurately measured 
Balmer decrement, can recover the intrinsic optical emission line 
luminosities to better than a factor of two in the majority of ordinary 
starforming and dusty galaxies. Our data shows that this remains 

12 Face-on galaxy orientations can also cause increased sca ttering into the 
line-of-sight and thus greyer dust curves i Rocha et all 20()8). 



true for galaxies with dust contents up to an Av of 4.4. However, 
while the use of a MW, /10.3 or Calzetti dust-curves for large sam- 
ples of galaxies will cause the least overall bias in the final dust- 
corrected emission line luminosities, the effect of a greyer dust- 
curve in individual galaxies should not be understated. A galaxy 
with a very typical measured Ha/H/3 ratio of 4.5 would have its 
Ha luminosity, and thus star-formation rate, corrected by a factor 
of around 2.5 using the average dust curve, whereas a factor of 7 
would be more appropriate if the majority of its dust were in the 
diffuse ISM. 



4.2 Fitting the dust curve 

Throughout the remainder of the paper we will focus on results us- 
ing the CF00 fi dust-curve as it allows greater flexibility than the 
traditional dust-curves in the calculation of errors, and the discus- 
sion of variations in dust-curves within the sample as presented in 
the previous section. Additionally, it is the only dust-curve specifi- 
cally motivated and justified for the correction of nebular emission 
lines. The MW dust-curve assumes a screen-like dust distribution, 
and the Calzetti dust-curve is intended for correcting the contin- 
uum light from galaxies. Rather than assuming n — 0.3 as derived 
by CF00, we can fit the [Ne] vs. Ha luminosities of the SSGSS- 
starforming galaxies to obtain the best-fit fi and intrinsic [Ne]/Ha 
ratio together with errors. We use a maximum-likelihood linear fit, 
accounting for errors on both quantities, and fix the intercept to be 
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Figure 7. The optical Balmer decrement vs. Silicate 9.7/rni absorption 
strength. The Balmer decrement is defined with an intrinsic ratio of B=2.86 
(3.1) for galaxies classed as starforming/composite (AGN) from their opti- 
cal emission lines. Thus, a Balmer decrement of 1 indicates no optical dust 
attenuation. The top axis indicates the corresponding optical ry, assum- 
ing the best-fit /jo.4 dust-curve derived in Section l4~2l The right-hand axis 
indicates the corresponding mid-IR 7y from the zero-age starburst mod- 
els of Rowan-Robinson & Efstathiou (2009). Overplotted as a dotted line 
is the best-fit relation to the SSGSS-starforming, composite and Balmer- 
strong galaxies as described in the text. Only those galaxies with Ha and 
H/3 detected at greater than 3cr and mean per-pixel-SNR> 2 between 5.4 
and 5.9/im are plotted. Symbols are: SSGSS starforming galaxies (green 
dots, 25/35), SSGSS AGN (green open triangles, 3/4), SSGSS composite 
galaxies (green open diamonds, 14/21), Balmer-strong galaxies (red dots, 
11/11), ULIRGs (purple stars, 6/6), Seyferts (blue triangles, 13/20). 

zcr Formal, statistical errors on the derived /i are at the level of 
5%. The errors obtained from bootstrap resampling with replace- 
ment of the data are more significant, and we quote these errors 
below, and use them to obtain errors on the dust correction formu- 
lae presented in Section [5Tl and l5"72l 

We find = 0.4±0.2, i.e. on average 40% of the optical depth 
at 5500A arises in the ISM of the galaxies. This is consistent within 
the errors with the fi — 0.3 derived by CF00 for their sample of 
local starburst galaxies. The intrinsic [Ne]/Ha ratio fit at the same 
time is g = 0.6 ± 0.1. The resulting curve is shown as the solid 
gray line in Figure|6]and given explicitly as: 



— = 4 1 = 0.6(A/5500)~ 

TV Ay 



L3 +0.4(A/5500)"°' 7 



(7) 



where A is in A. Given a measured Ha/H/3 flux ratio, the optical 
depth at any wavelength can be obtained from: 



3.22 In 



2.86 



(8) 



J(Ha)/f(H/3)J r v 

where t\/tv is given in Eqn.lT^A = 1.086ta is the attenuation in 
magnitudes, and a Case B recombination ratio of 2.86 is assumed. 

4.3 Dense vs. diffuse dust 

Extinction in the mid-IR continuum is dominated by the stretching 
and bending modes of amorphous silicate grains with two strong 
features at 9.7/im and 18/im. Importantly for this section, it is evi- 
dent from the many theoretical models that any significant silicate 
absorption, as found in the ULIRGs and many of the Balmer-strong 

13 Allowing the intercept to be non-zero results in the same fitted parame- 
ters within the errors. 



galaxi es, implies optical extinction at levels of tv ~ 50 — 100 
(e.g. Isirocky et alj|2008l : iRowan-Robinson & Efstathioull2009h . It 
is therefore impossible for ionising UV continuum and associated 
line radiation to pass through these clouds directly. 

In Figure|7]we plot the Balmer decrement [Ha/(B x H/3) where 
B is the unattenuated ratio] against the 9.7/xm Silicate absorption 
strength for all four samples. For H II regions we set B = 2.86 
for Case B reco mbination at lO.OOOK and e lectron density of 
10 2 - 10 4 /cm~ 3 (IOsterbrock&Ferlandll2006T) . For AGN we set 
B — 3.1 JtCewlev et al.ll20061) . Surprisingly, we see a correlation 
between these two measures of dust content for most galaxies in 
our samples. However, it is clear that they do not trace the same 
dust because the tv measured from the Balmer emission lines (up- 
per axis) is an order of magnitude lower than that indicated by the 
mid-IR silicate absorption strengths (right-hand axis, from zero-age 
starburst models of Rowan-Robinson & Efstathiou 2009). The dot- 
ted line is the best-fit squared relation to the SSGSS-starforming, 
composite and Balmer-strong galaxies: 



si 



9.7/jm 



= -0.17 



/(Ha) 



(9) 



In the case of starforming galaxies without a strong AGN, 
the mid-IR continuum predominantly arises from dust grains sit- 
uated within the H II regions re-radiating the light from young and 
hot stars. It is this dust continuum that is absorbed by the Silicate 
grains along the line-of-sight between the light source and the ob- 
server. Balmer emission lines also come from ionised gas within 
Hll regions. Thus, the observed correlation between the Silicate 
absorption strength and Balmer decrement implies either: (i) a cor- 
relation between the amount of dust in dense stellar birth clouds 
which cause the Silicate absorption, and the amount of dust in more 
diffuse stellar birth clouds which cause the absorption of the opti- 
cal emission lines; or (ii) a self-similar structure for H II regions in 
which the Balmer lines that reach the observer have passed through 
a fixed fraction of the dense dust clouds. 

There are two notable exceptions where the optical Balmer 
decrement fails to predict the presence of an extremely signifi- 
cant attenuation in the mid-IR: the ULIRGs 15250 and 08572. 
iNardini et all {2009) find that both of these ULIRGs have signif- 
icant contributions from AGN to their bolometric luminosity (53 
and 87% respectively), using a spectral decomposition between 5 
and 8^tm. Models have shown that extremely strong silicate ab- 
sorption, as observed in these two ULIRGs, requires a steep tem- 
perature gradient within the obscuring cloud, which can only be 
obtained where the light source responsible for heating the dust 
that causes the mid-IR continuum is deeply buried in a smooth 
thick shell of dust th at is both geometrically and optically thick 
jLevenson et alj2 007). This is again consistent with the source be- 
ing an AGN, although a nuclear starburst may suffice. A clumpy 
distribution of dust le ads to much shallower absorption or emission 
(Nen kovaetail 12002). Our emission line analysis is also consis- 
tent with these galaxies harbouring a buried AGN, they have an 
increased [Ne]/Ha ratio and lie in the composite region of the BPT 
diagram. While the position of these galaxies on the BPT diagram 
is not in itself conclusive proof that these two ULIRGs harbour 
an AGN, if their position is caused by the AGN then it indicates 
that ionising radiation from the AGN escapes to ionise the NLR, 
despite the dense dust cloud that presumably surrounds the central 
continuum source. Two dusty Seyfert galaxies and one ULIRG with 
a Seyfert nucleus (Mrk273) also lie above the relation defined by 
the SSGSS and Balmer-strong galaxies, further supporting the idea 
that dust along the line-of-sight to a deeply buried AGN can lead 
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to stronger silicate absorption, relative to nebular line absorption, 
than normal H II regions. 

We are not aware of any models which simultaneously pre- 
dict the strength of silicate absorption and attenuation of the opti- 
cal emission lines. Qualitatively Figure|7]suggests a two component 
model is required to describe the dust distribution in galaxies/ AGN. 
The lower branch is consistent with silicate absorption arising in the 
clumpy distribution of stellar birthclouds, with the densest clouds 
causing some silicate absorption, and the Balmer emission lines 
originating from the same star formation regions but further out in 
the clouds. The upwards scatter of some sources from this branch 
is consistent with the dust being concentrated in a dense shell-like 
distribution obscuring a single nuclear source, most likely an AGN 
in the cases studied here, but also conceivably a strong nuclear star- 
burst. 



4.4 Galaxies with an AGN 

As explained in Section 12.11 a dust curve and accurate Balmer 
decrement that adequately corrects H II region emission lines for 
dust attenuation, will not necessarily do the same for lines from 
the NLR of AGN. Particularly in composite objects, a higher frac- 
tion of Balmer lines will arise from within the H II regions of the 
galaxy compared to lines such as [OIII] which are stronger in NLRs. 
Therefore a form of differential extinction might arise between 
the Balmer and higher ionisation lines, and a corresponding bias 
in the correction of AGN NLR lines for dust using the measured 
Balmer decrement. Additionally, the form of the dust curve may 
not be the same for light originating from Hll regions and from 
an AGN NLR, due to the different relative geometry between the 
light source and the dust cloud, and the different balance between 
dense and diffuse dust attenuating the lines. The mid-IR [OIV] line 
at 25.8 9 nm has been shown to be an accurate indicator of AGN 
power jMelendez et al]2008h , and the [OlII]/[O lV] ratio a possible 
measure of AGN attenuation JHaas et alj|2005h . In this section we 
test the accuracy of dust-corrected L[OlII] as a measure of AGN 
luminosity in our sources by comparing to the mid-IR [OIV] line. 

In Figure [8] we plot dust-uncorrected and dust-corrected 
L[OlII] vs. L([OlV] A26/im) for the Seyferts, the 3 Balmer-strong 
galaxies and the 2 ULIRGs with measured [OIV] lines. Note that 
the correlation coefficient is similar for the Seyferts whether the 
dust correction is performed or not. These Seyferts have very small 
Balmer decrements so the dust correction has only a small effect. 
Because some of these objects are classified as composite AGN- 
starformation by their optical line ratios, and [OIII] is also pro- 
duced by starformation, in the right ha nd panel we apply a cor- 
rection for this following the method of IWild et al.l bOlCf) . based 
on the expected [OlII]/Ha ratio for metal-rich starforming galax- 
1C [OIV] is also produced in H II regions, although the contam- 
ination is expected to be much smaller than for [OIII] due to the 
higher ionisation level. Clearly, there is a tight correlation between 
dust-corrected L[QII I] and L[OlV] for the Seyfert sample (see also 
LaMassa et al. 2010). All of the dusty Balmer-strong galaxies and 
ULIRGs with measured [OIV] fall within a factor of 3 of the re- 
lation defined by the Seyferts. Only the upper limit for Arp220 
is marginally inconsistent with the relation defined by the Seyfert 
galaxies. Taking the remainder of the Balmer-strong galaxies and 
stacking the SH spectra in luminosity units shows a clear [OIV] 

14 We note that such a correction would not be appropriate for the ULIRGs 
if shocks were a significant contributor to the total Ha luminosity. 



line, which we fit together with the neighbouring Fell line using 
a double Gaussian. The corresponding stacked L[OlV]-L[OlII] lu- 
minosity also lies exactly on the Seyfert relation, once contamina- 
tion of L[OlII] from star formation has been accounted for. Errors 
are estimated from the standard deviation of the luminosity on re- 
moving one spectrum at a time from the stack. This strong corre- 
lation between [OIV] and [OIII] dust corrected using the observed 
Balmer decrement and dust-curve appropriate for H II regions, for 
even dusty composite galaxies, shows that there is no evidence for 
significant differential extinction between lines emitted from the 
NLR and H II regions of a galaxy. Although slightly surprising, this 
may imply that the majority of dust attenuation suffered by NLR 
emission lines comes from the diffuse dust and H II regions in the 
intervening galaxy, through which they must pass on their way to 
the observer. 

Unfortunately neither of the two ULIRGs with extreme sili- 
cate absorption (08572 and 15250) have measured L[OlV], so we 
cannot test directly whether optical and mid-IR NLR lines agree 
even in the presence of extreme, dense obscuration of the central 
light source. Such a study may be possible with the large samples 
of ULIRGs now available in archives. 



5 DISCUSSION 

The potential of dust attenuation and complicated dust geometries 
to render useless optical observations of extreme galaxies is some- 
times used to justify rest-frame IR observations. In the previous 
sections we have used the correlation between the unattenuated 
mid-IR and attenuated optical nebular emission lines to argue that 
a single dust curve can accurately correct optical emission lines for 
dust attenuation in a large range of different galaxies. The correc- 
tion is accurate to within a factor of 2 in most cases, and to at least 
a factor of 4 for even the most heavily obscured ULIRGs. 

Here, we discuss the implications of our work for deriving 
SFRs and black hole accretion rates from the optical emission lines 
Ha, [On], and [OIII]. We present useful relations for dust correc- 
tions at the wavelengths of each of these lines. 

5.1 Star formation rates 

The left and right panels of Fig. [9] show the fractional attenuation 
(in flux units) as a function of observed Balmer decrement at Ha 
and [Oil], for each of the different dust curves discussed in Section 
|4] Using the best-fit /iq.4 dust curve, the recommended correction 
at Ha is: 

L[Ha] corr = L[Ha] obs ( ( Ua /^U S \ ' " ' (1Q) 
and at [Oll£l 

L[OII] corr = L[OII] obs ^ ( HQ /H/?)o bs ^ ■ " - (U) 

where B is the unattenuated ratio as discussed above and errors are 
propagated from the errors on the best-fit fj,, While for H II galaxies 
and strong Seyferts the choice of B is clear, for composite galaxies 
the unknown B may cause small systematic uncertainties. However 

15 Note that this is an extrapolation of the dust-curve beyond the wave- 
length region directly tested in this paper. 
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Figure 8. [OIII] A5008A vs. [OIV] A26/im luminosities (erg/s) for the Seyfert sample (filled blue triangles, 20/20), ULIRGs (purple stars, 2/6), and the Balmer- 
strong galaxies (red circles, 3/11) with [OIV], [OIII], Ha and H/3 detected at > 3<r. The stack of the remaining Balmer-strong galaxies with non-detections 
in [OIV] is shown as a red open circle. Left: L[OIII] is not corrected for dust attenuation, or for contamination due to star formation. Note the different 
x-axis range in this panel. Center: L[OIII] is corrected for dust attenuation but not for contamination due to star formation. Right: L[OIII] of the ULIRGs 
and Balmer-strong galaxies has been corrected for contamination from star formation using the method described in Wild et al. (2010). Dust attenuation is 
calculated using the best-fit dust-curve derived in Section |4~21 Overplotted as a dashed line is the best-fit linear relation to the Seyfert sample only, the dotted 
lines indicate a factor of three above and below this relation. For clarity, error bars are omitted in the right hand panel. 




Figure 9. Intrinsic to observed flux ratios as a function of observed Hct/H/3 line ratio for different dust curves, at wavelengths corresponding to 
[Oi l] AA3 727, 3730(left), [OIII] A5008A (center) and Ha (right). The lines are (as given in Figure|6}: the SMC extinction curve with an R v = 3.1 (dotted 
Pell 19921): a MW extinct ion curve with Ry = 3.1 (dashed black, IO'Donnelj |i994): a different MW ex tinction curve wi th Rv =3.2 (dashed grev. lSeatonl 
1979l ;lNandvetalJll975t) the /in. 3 and /in. 4 attenua tion curve fo r emission lines (solid black and grey, [Wild et al.f 2007); the Calzetti attenuation curve for 
starburst continua with Ry = 4.05 (dash-dot ICalzetti et alJ ioOO): the A~ 0,7 attenuation curve for starburst continua from CF00 (dash-triple-dot). Note the 
changing scale on the y-axis. 



as Hll regions have relatively stronger Ha emission, a value of 
B = 2.86 is probably close to reality. 

To derive a star formation rate from optical nebular emission 
lines, two steps are required. This paper has focused on the first 
step: to correct the line for dust attenuation. The second step is to 
convert the dust-corrected line luminosity into a star formation rate. 
While our results show that dust-corrected Ha luminosity traces 
[Ne] luminosity, this is not enough to confirm the reliability of Ha 
as an indicator of star formation rate in extremely dusty objects. 
This is because dense dust clouds can have high optical depths even 
in the mid-IR, as indicated by the strong Silicate absorption in some 
galaxies, and thus even the neon lines may be attenuated in some 
objects. In addition, the most extreme compact, dusty galaxies may 
even have a non-negligible fraction of their ionizing photons ab- 
sorbed by dust and not gas (and thus not traced by Ha). Encour- 
agingly, Ho & Keto (2007) have already shown how [Ne] lumi- 



nosity correlates strongly with both the total infrared luminosity 
and Brackett-a luminosity in star-forming galaxies. Their galaxy 
sample includes a handful of ULIRGs, showing that these relations 
hold even at extreme luminosities. This result provides circumstan- 
tial evidence that the neon luminosity is tracing star-formation even 
in the dustiest of galaxies. In contrast, Farrah et al. (2007) find [Ne] 
luminosity in ULIRGs to be deficient by 0.4dex from the Ho & 
Keto relation, which they ascribe to extinction of the [Ne] lines. In 
a companion paper we plan to address this issue directly by com- 
paring the star formation rates derived from the full SEDs, with that 
derived from the emission lines. 



5.2 Black hole accretion rates 

As with star formation rate, the derivation of black hole accretion 
rates (BHAR) from nebular emission line luminosities requires two 
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steps: correction for dust emission discussed in this paper, and then 
conversion into BHAR via a bolometric correction. 

For the best-fit /Uo.4 dust curve, the relevant correction for dust 
attenuation at [OIII] is: 



L[OIII], 



L[OIII] 



ga/Hgg 
B 



e + 0.40 
-0.32 



(12) 



where B is the unattenuated ratio as discussed above. The correc- 
tion caused by the different dust curves is presented in the central 
panel of Figure|9] 

The exponent for the /10.4 dust curve is larger than that used 
in some AGN literature and can lead to significant differences 
in estimated bolometric AGN luminosity and accretion ra tes as 
a fraction of LEdd (see also the discussion in Netzer 2009). Us- 
ing a smaller exponent than appropriate will cause an underesti- 
mate in L[OlII] C orr which will vary with dust content. For exam- 
ple, for the relatively dust-free sample of Seyfert galaxies stud- 
ied in this work, the difference between the fio.4 dust curve ex- 
ponent of 3.55, and 2.9 4 (from a MW extincti on law) used in e.g. 
lBassanietai]jl999h and lLamastra et al" I J2009T) leads to differences 
ranging from 5 to 30%. The potential increase in the under- (or 
over-)estimation of the line luminosity with increasing dust con- 
tent because of an inappropriate dust-curve has important impli- 
cations for understanding whether the bolometric correction for 
L[OlII] changes systematically with L[0III], due to the correla- 
tions between SF R of the host galaxy and the luminosity of the 
central AGN (e.g . i Kauffmann etai1l2003l : ISchweitzer etal]|2006l : 
IZheng et al.ll2009l:lNetzeill2009h and the SFR and dust content of 
the host 1 da Cunha et al. 2010h . Our results do not explicitly con- 
strain which dust-curve is most appropriate for Seyferts, however, 
we would like to emphasise the importance of including the errors 
caused by the dust-attenuation correction when studying emission 
line trends. 



5.3 AGN detection 

Even in the two ULIRGs in our sample where there is evidence for 
a heavily buried central nucleus, we would like to emphasise that 
the optical emission line ratios still indicate the possible presence 
of an AGN, albeit not conclusively. This shows that even where the 
AGN is deeply buried, some radiation may still escape to ionise 
the NLR clouds. This has important implications for using optical 
emission line ratios to identify AGN in objects such as ULIRGs, 
and will become increasingly important at high-z. Unfortunately 
neither of the two ULIRGs with strong silicate absorption in our 
sample have measured [OIV], and therefore we cannot directly ver- 
ify that their position on the BPT diagram, and therefore their cor- 
rection of [OIII] for contamination from star formation, is correct. 
The specific study of a larger sample of ULIRGs to compare the 
bolometric AGN luminosities derived from mid-IR lines, mid-IR 
continuum and dust- and star-formation corrected high-quality op- 
tical emission lines would now be possible with samples in the 
SDSS and Spitzer archives. 



6 CONCLUSIONS 

The purpose of this paper was to verify, or otherwise, the accu- 
racy of fundamental quantities derived from optical emission lines, 
particularly in the case of significant dust attenuation of the lines. 
In particular, optical emission lines provide powerful measures of 
star formation rate and black hole accretion rate. We compared two 



"dusty" galaxy samples, including 6 ULIRGs, to a sample of "ordi- 
nary" galaxies from the SSGSS all with high quality SDSS optical 
spectra and mid-IR spectra from the Spitzer Space Telescope IRS 
spectrograph. The results are extremely encouraging for using rest- 
frame optical spectra to study dusty galaxies with an Av < 4.5. 
Our results are promising even for galaxies as dusty as ULIRGs, 
although a much larger sample should be studied covering the full 
range of this diverse population, in order to verify our conclusions. 
Our results have important implications for future high-redshift 
spectroscopic surveys, where galaxies have higher star formation 
rates and dust contents may be correspondingly higher. 

We have compared several different dust curves used in the 
literature to correct optical emission lines for dust attenuation. Our 
results favour a dust curve which lies close to the MW extinction 
curve, indicating that a significant fraction of the dust which atten- 
uates nebular emission lines has a screen-like geometry. This is in 
ag reement wit h the model for emission li n e atte nuation presented 
by IWild et all J2007t) and Ida Cunha et al.l (2008) and based upon 
the 2-component CF00 dust model. In this model, the screen-like 
extinction arises from the dense stellar birth clouds which surround 
the hottest stars responsible for ionising the gas. We find a best- 
fit dust-curve with on average 40% of the optical depth at 5500A 
arising from diffuse ISM dust. The /10.4 dust curve is not signif- 
icantly different from the Calzetti dust-curve in the optical wave- 
length regime, however, we focus our results and analysis on the 
former which was intended for emission lines and allows greater 
flexibility, rather than the latter which was derived from galaxy con- 
tinua. 

Particular conclusions with regard to measuring star formation 
from optical emission lines in the presence of dust: 

• The MW, Calzetti, /io.3 and best-fit /io.4 dust curves all allow 
accurate correction of Ha (to better than a factor of 2), and there- 
fore accurate calculation of star formation rates, for galaxies with 
observed Ha/H/3 ratios of as much as 10 (ry < 4.5). 

• For only 2/6 ULIRGs with Ty < 6 does the optical Ha differ 
from the mid-IR [Ne] lines by more than a factor of 2, but then 
only by a factor of 4. In these two cases (08572 and 15250) we 
have argued that the observed excess of [Ne] is in fact caused by 
a buried AGN, rather than u ncertainties in the d ust correction (in 
agreement with the results of lNardini et alj|2009l) . 

• The majority of galaxies in our samples favour a dust curve 
close to the MW screen-like dust extinction curve. This suggests 
that much of the extra dust in the unusually dusty objects is located 
in dense clouds, rather than in the diffuse ISM. 

• There is tentative evidence for a range in the shape of the dust 
curves in ordinary starforming galaxies. For the CF00 dust model, 
this range translates into a variation in the fraction of optical depth 
arising in the diffuse ISM of between around 30% and 90%, with 
the remaining fraction arising in dense birth clouds. However, the 
variation in the intrinsic [Ne]/Ha ratio with metallicity, the issue 
of relative aperture correction between the IRS and SDSS data, and 
the large statistical errors on the H/3 emission line, prevent us from 
quantifying this scatter more precisely. 

Particular conclusions with regard to identifying AGN and 
measuring their accretion rates from optical emission lines in the 
presence of dust: 

• The best-fit /10.4 dust curve and measured Balmer decrement 
allows accurate correction of [OIII] for dust attenuation, to within 
a factor of 3. This is slightly surprising, and suggests that the ef- 
fects of variations in the relative geometry of the source and dust, 
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and differential extinction between the Balmer and higher ionisa- 
tion lines, are minimal. 

• The method to separate the contributions to optic al emission 
lines f rom both star formation and AGN presented in IWild et al.l 
(2010) tightens the relation between L[OlV] and L[0III] for com- 
posite objects. This supports the use of the position of galaxies on 
the optical BPT line ratio diagram to measure the amount of L[OlII] 
originating from the AGN, and therefore the estimates of black hole 
accretion rates from L[0III] corrected for dust attenuation and star 
formation contamination. 

• All ULIRGs in our sample are classified as AGN or compos- 
ite objects from their optical emission line ratios. While the com- 
posite line ratios may be caused by shocks, we have argued that 
other observational results favour the presence of a buried AGN 
contributing significantly to the mid-IR continuum. It is therefore 
possible that sufficient ionising flux escapes from the central source 
to ionise the NLR, even in the case of heavily buried nuclei such 
as in ULIRGS 08572 and 15250. Unfortunately these two ULIRGs 
do not have measured L[OlV] preventing us from verifying our cor- 
rections to L[0III] for dust attenuation and star formation contami- 
nation in such extreme cases. In the case of ULIRGs, whose prop- 
erties vary widely, clearly a much larger sample should be studied 
before firm conclusions can be drawn. 
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APPENDIX A: SAMPLE INFORMATION 

APPENDIX B: Balmer-strong GALAXIES 

While three of the samples studied in this paper have been pre- 
sented elsewhere, the Balmer-strong sample has not. Here we sum- 
marise its selection criteria and the mid-IR observations. 

Bl Selection 

The galaxies were selected to have: 

• Observed Ha/H/3 flux ratios greater than 8 

• Strong Balmer absorption lines compa red to their 4000 A 
break strength [PC2 > 0.05 x PCI + 0.4. IWildet all J2007h l. 
Although they were not selected on a single Balmer absorption line 
equivalent width, which have considerably larger errors than the 
spectral indices used here, most have US a > 5A and D n (4000) < 
1.4. 

• Accurately measured emission lines: a signal-to-noise > 
3a for all four BPT lines, together with a H/3 flux > 40 x 
10" 17 erg/s/cm 2 . 

• A contribution to their emission lines by an AGN, i.e. they lie 
above the empirically derived demarcation line of Kauffman n et al.l 
d2003h on the BPT diagram. 

• z < 0.07 so that the SDSS spectrum, taken through a 3"fibre, 
samples only light coming from the central few kpc of the galaxy. 

B2 Mid-IR observations 

IRS observations in both high and low resolution were carried out 
using staring mode during program 40330 (PI Timothy Heckman). 

We used SPICE to extract the ID spectra from the low- 
resolution basic calibrated data (BCD) pipeline products, using the 
point-source and optimal extraction options. The two ID spectra 
(i.e. one for each nod) were combined using a weighted mean. For 
the high-resolution spectra, we first identified and interpolated over 
additional rogue pixels using a custom IDL routine and the rogue 
pixel masks suitable for our campaign numbers (41-46). This pro- 
cedure was repeated for the dedicated off-source background ob- 
servations. All 3 science exposures were combined using a sim- 
ple arithmetic mean, resulting in two 2D spectra. These were sky- 
subtracted and the ID spectra were extracted using SPICE with 
the point-source and optimal extraction options. For each nod, the 
edges of each order were inspected and trimmed and the orders 
combined using a weighted mean where wavelengths were covered 
by more than one order. The ID spectra of each nod were subse- 
quently combined using a weighted mean. 
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Name 


RA (deg) 


Dec (deg) 


specobjid 


rowindex 


z 


PSB01 


181.225 


52.395 


7484858593665351 68 


284940 


0.0632 


PSB02 


169.600 


56.036 


256085720727289856 


296570 


0.0680 


PSB03 


173.251 


60.275 


268189274979958784 


315471 


0.0644 


PSB04 


167.291 


45.857 


404988882665340928 


491844 


0.0636 


PSB05 


136.004 


1.458 


132797726868897792 


105860 


0.0535 


PSB06 


120.199 


37.729 


213582458531086336 


232315 


0.0416 


PSB07 


339.726 


13.194 


208235580251701248 


223368 


0.0628 


PSB08 


128.311 


32.583 


262560681922920448 


306604 


0.0558 


PSB09 


239.431 


27.465 


391758540508758016 


472114 


0.0316 


PSB10 


139.864 


33.791 


448335754997792768 


536057 


0.0193 


PSB11 


158.162 


12.178 


450306505485320192 


539391 


0.0329 



Table A2. The ULIRG sample, where SDSS spectra exist. Columns are: name as given in Armus et al, (2007) [1], RA [2], Dec [3], SDSS specobjid [4], 
SDSS-MPA rowindex [5], redshift from SDSS [6]. 



Name 


RA (deg) 


Dec (deg) 


specobjid 


rowindex 


z 


IRAS 08572+3915 


135.106 


39.065 


337714854909444096 


394321 


0.0579 


IRAS 12112+0305 


183.441 


2.811 


146028589375029248 


126942 


0.0730 


IRAS 15250+3609 


231.748 


35.977 


394574695615168512 


476759 


0.0553 


Arp220 


233.738 


23.504 


609061542893191168 


755340 


0.0181 


Mrk273 


206.176 


55.887 


372055066565148672 


439991 


0.0373 


UGC5101 


143.965 


61.353 


137019791402598400 


111495 


0.0393 



Table A3. The Seyfert sample. Columns are: name as given in lLaMassa et alj l2010l) ril. RA [2], Dec [3], SDSS specobjid [4], SDSS-MPA rowindex [5], 
redshift from SDSS [6]. 



Name 


RA (deg) 


Dec (deg) 


specobjid 


rowindex 


z 


2MASXJ08244333+2959238 


126.180 


29.990 


339966523504328704 


397741 


0.0250 


CGCG064-017 


149.812 


12.988 


491120230488080384 


603606 


0.0340 


Mrk0609 


51.356 


-6.144 


129701503752470528 


101843 


0.0340 


NCG0291 


13.375 


-8.768 


185153160128495616 


190581 


0.0190 


NGC5695 


219.342 


36.568 


389226545811554304 


467780 


0.0140 


2MASXJ 1 5522564+2753435 


238.107 


27.895 


391758541012074496 


472216 


0.0740 


2MASXJ12384342+0927362 


189.681 


9.460 


347285137836212224 


409636 


0.0830 


2MASXJ08035923+2345201 


120.997 


23.756 


356292212926971904 


415287 


0.0290 


2MASXJ 1 6 1 64729+37 1 6209 


244.197 


37.273 


297464197648744448 


366718 


0.1520 


2MASXJ1 1570483+5249036 


179.270 


52.818 


248485858699640832 


284823 


0.0360 


2MASXJ 12 1 83945+4706275 


184.664 


47.108 


408648327781941248 


498207 


0.0940 


SBS1 133+572 


173.954 


56.952 


369240319515951104 


435171 


0.0510 


SBS1609+490 


242.716 


48.911 


175301007746531328 


174339 


0.0450 


IC0486 


120.087 


26.614 


261716087635181568 


305068 


0.0270 


CGCG242-028 


170.755 


47.052 


405833433537839104 


493622 


0.0250 


Mrkl457 


176.840 


52.450 


247922896524017664 


283914 


0.0490 


2MASXJ 1 34632 1 7+6423 247 


206.634 


64.390 


170234565783191552 


166011 


0.0240 


CGCG2 18-007 


200.952 


43.301 


387537584101785600 


464960 


0.0270 


2MASXJ10181928+3722419 


154.580 


37.378 


401892408733728768 


486221 


0.0490 


2MASXJ1 1 1 10693+0228477 


167.779 


2.480 


143495710086529024 


122473 


0.0350 



© 0000 RAS, MNRAS 000, 000-000 



